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We present observations of the H91q; recombination line emission towards a sam- 
ple of nine Hll regions associated with 6.7-GHz methanol masers, and report arcsecond- 
scale emission around compact cores. We derive physical parameters for our sources, 
and find that although simple hydrostatic models of region evolution reproduce the 
observed region sizes, they significantly underestimate emission measures. We argue 
that these findings are consistent with young source ages in our sample, and can be 
explained by existence of density gradients in the ionised gas. 
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1 INTRODUCTION 

Ultracompact (UC) Hll regions are pockets of ionised hy- 
drogen that form around massive stars in the earliest stages 
of their evolution. Together with massive bipolar outflows, 
strong far infrared emission by dust, and the presence of 
molecular masers, they are indicative of massive star forma- 
tion in its earliest stages. The simple model of Hll region 
evolution jSpitzer 1978^ does not explain many of their ob- 
served properties, such as the frequent occurrence of non- 
spherical morphologies and the lifetime paradox. This life- 
time problem, first noted by Wood & Churchwell ( [158.91 , is 
that the number of observed UC Hll regions exceeds that 
predicted from their dynamical expansion timescales by two 
orders of magnitude, given the accepted massive star forma- 
tion rate in the Galaxy. A number of modifications and en- 
hancements to the basic model have been suggested, includ- 
ing the work of Dyson et al. ( 1995 1 , HoUenbach et al. 1199411 , 
Tenorio-Tagle 119791 and van Buren et al. II199U1 . Franco 
et al. inMIl, Arthur fc Liza no llTWIl . Keto ll^IMl. Thether- 
mal IjDe Free et al. 1995at and turbulent IXie et al. 199611 
pressure confinement models are appealing due to their de- 
pendence on the ambient conditions observed to commonly 
exist in molecular clouds. 

De Free et al. Ill995all suggested thermal pressure con- 
finement as an explanation of the lifetime paradox. They 
noted that when Wood & Churchwell proposed the lifetime 
problem in 1989, the molecular medium surrounding the 
UC Hii regions was thought to have temperatures ~ 25 K 
and densities ~ 10® cm~^. More recent observations indicate 
T~ 100 K and n~ lO'^ cm"^. The resuhing 400 x increase 



in thermal pressure limits the expansion of the Stromgren 
sphere. A weakness of this model, noted by Xie et al. Ill996ll . 
is the exceedingly high emission measures that it predicts 
(~ 10^° pc cm~®), which are more than two orders of mag- 
nitude greater than the values typically observed. 

Hierarchical density and temperature structures are 
known to exist within star- forming regions, with hot cores 
embedded in larger, less dense molecular clumps which 
themselves are within still larger and less dense molecu- 
lar clouds. The densities decrease by approximately an or- 
der of magnitude in going from core to clump and again 
from clump to cloud (Ces aroni et al. 1994t . In a seminal 
work. Franco et al. 119901 showed that these density in- 
homogeneities are important for Hll region evolution. That 
the hierarchical structure of molecular clouds plays an im- 
portant role in Hll evolution is supported by the marked 
similarity in the ionised and neutral gas density struc- 
tures observed within molecular clouds (IKim fc Koo 199'6l 
IKim fc Koo 20021 IKoo fc Kim 20031 . 

The large thermal molecular line and recombination line 
widths observed towards many UC Hll regions suggest that 
significant turbulent motions are present, probably of mag- 
netic origin | |Garci'a-Segura fc Franco 1996| l. This led Xie 
et al. 119961 to suggest that turbulent pressure is the domi- 
nant mechanism to restrict the expansion of an Hll region. In 
contrast to thermal pressure confinement, the assumed den- 
sities are lower, resulting also in lower emission measures. 
During much of the expansion phase, however, the turbu- 
lent pressure is expected to play a lesser role than thermal 
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pressure in confining the Hll region. This is discussed briefly 
in section Pi. 1.1 1 

Icke (1979) investigated the formation of Hll regions in 
non-homogeneous media and was able to explain some non- 
spherical morphologies. Observational evidence obtained 
in the last decade, however, suggests that many Hll re- 
gions have compact cores within diffuse, arcminute-scale 
extended emission IIKurtz et al. 19991 IKim fc Koo 200111 . 
Other sources exhibit this to a smaller degree — the so- 
called core-halo morphology; see Wood & Churchwell II1989II 
and Kurtz et al. 1199411 . A study of the compact and ex- 
tended radio continuum emission and radio recombination 
lines (RRLs) from eight Hll regions known to be associated 
with 6.7-GHz methanol masers has been undertaken. The 
RRL analysis is presented here, while details of the con- 
tinuum observations can be found in Ellingsen, Shabala & 
Kurtz (hereafter ESK05). In sectionHwe briefly out- 

line our observations. The results are presented in section |3 
and these are compared with a simple model in section 2] A 
discussion of our findings is presented in section |3 



2 OBSERVATIONS 

Eight UC Hll regions with associated methanol maser emis- 
sion were observed with the Australia Telescope Compact 
Array (ATCA) on 1999 July 10 and 11. Both continuum 
and recombination line emission were observed in the ATCA 
75QD array, with an angular resolution of 7" and largest de- 
tectable angular scale of ~ 50". Details of the continuum 
observations are found in ESKQ5. 

The H91a recombination line {uq = 8.58482 GHz) 
was observed with an 8-MHz bandwidth and 512 spec- 
tral channels, giving a frequency resolution of 15.625 kHz 
(0.529 km s"'^), and total velocity coverage of 270 km s~^. 
The data were further smoothed by frequency averaging over 
four or eight channels. All sources were observed together 
with associated secondary calibrators immediately before 
and after each on-source observation. The primary calibra- 
tor PKS B1934-638 was observed each day to calibrate the 
flux density scale. 

The observations were made with the array in the 750D 
configuration, with a minimum baseline of 31 m and max- 
imum of 719 m. Because the primary aim of these obser- 
vations was to look for extended emission associated with 
the sources, baselines to the 6 km antenna were not used 
in order to maximize sensitivity to large-scale structure. A 
summary of the fields observed is given in Table 1. 



3 RESULTS 

H91a emission was detected in six sources. The observed 
line profiles and Gaussian fits obtained after continuum sub- 
traction are shown in figure As we are only interested in 
emission around the RRL peak, non-zero baselines were used 
in some sources for a better fit in those regions. The flux 
densities per channel were obtained by integrating over the 
area of the UC Hll region (typically around 5" x 5"; see Ta- 
ble|5J. No H91a emission was detected from G336.40— 0.25, 
G339.88-1.26 or G345.0 1-1-1. 79. These three regions have 



low continuum brightness so the expected LTE line bright- 
ness is at or below the image noise level. 

The fit parameters are given in Table |5| These were 
achieved using a standard non-linear least squares fitting 
routine, and resulted in unrealistically low error estimates 
for the fit parameters. In order to determine a realistic 
uncertainty, we made Monte-Carlo simulations, using the 
same routine to fit Gaussian profiles with known parame- 
ters, plus white noise of various amplitudes. Comparison of 
the methanol maser peak velocities from Table Q with the 
H91a velocities from Table|5|shows approximate agreement 
(indicating an association between the masers and the star 
formation region) but sufficient difference in some cases to 
suggest that the masers may not be directly linked to the 
ionised gas. The detection is only marginal for G309. 92-1-0. 48 
and thus no formal parameters were derived for this source. 

3.1 Arcsecond-scale emission 

3.1.1 Moment maps 

Image cubes for the six sources that exhibited H91a emis- 
sion were analysed for the presence of arcsecond-scale ex- 
tended emission near the UC Hll regions. Two complemen- 
tary methods were employed for this analysis. 

The first of these involved examining cross-cuts through 
the image cubes in various position-velocity planes. Posi- 
tional cross-cuts were taken to run through the peak bright- 
ness position at a range of angles. Typically, the cuts were 
made at constant right ascension or declination or along 
a line joining the compact region to features of poten- 
tial interest. Part a (left plots) in Figures |21 - |S] show the 
position-velocity plots of integrated flux density for appro- 
priate cross-cuts in sources exhibiting significant arcsecond- 
scale extended emission. Position angle is calculated anti- 
clockwise from a cut in constant declination. Position is de- 
fined as the offset from the plane passing through the point 
of maximum H91a emission. If both compact and extended 
components are present and physically associated, provided 
no strong shocks are present, their systemic velocities should 
be approximately equal or change smoothly between the two 
positions. Therefore, taking a cut along the line joining them 
would result in emission at the same peak velocity, but offset 
positionally by the distance between the two components. 

Another way of determining whether any two compo- 
nents are likely to be associated is by plotting the first mo- 
ment (a fiux-weighted velocity mean across the cube) distri- 
bution across the region of H91a emission. This is shown in 
greyscale in part b (right plots) of each figure. Superimposed 
are the continuum contours of observations also made in the 
750 D array configuration. 

Clipping levels used in Fig. |21 - HI were set to 2a, cor- 
responding to 6 - 11 mjy beam~^ depending on the source 
(see Table . Lower density gas outside the UC Hll cores 
is indicated in the first moment plots by presence of more 
emission at velocities close to the peak H91a velocity of 
the source (given in Table|2l. In position- velocity plots, this 
corresponds to "peaks" and/or "troughs" (depending on the 
location of extended emission) in position at the peak RRL 
velocity. Evidence of such emission is seen in G 318.95—0.20 
(FigurelU, G 328.814-0.63 (FigureOJ and NGC 6334 E (Fig- 
ure [IJ and F (Figure 01 components. Less arcsecond-scale 
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Source Methanol peak Right Ascension DecHnation Associated Central Wlsr Avg noise per spectral channel 

(kms-i) (J2000) (J2000) IRAS Source (km s"!) (mjy beam-i) 



G 308.92+0.12 


-54.5 " 


13:43:02 


-62:08:51 


13395-6153 


21.8 


2.3 


G 309.92+0.48 


-59.6 


13:50:42 


-61:35:10 


13471-6120 


21.6 


2.0 


G318.95-0.20 


-34.7 


15:00:55 


-58:58:42 


14567-5846 


18.3 


3.1 


G 328.81+0.63 


-44.0 = 


15:55:48 


-52:43:07 


15520-5234 


14.7 


3.9 


G 336.40-0.25 


-85.3 " 


16:34:11 


-48:06:26 


none 


11.2 


3.3 


G 339.88-1.26 


-38.7 


16:52:05 


-46:08:34 


16484-4603 


9.3 


2.8 


G 345.01+1.79 


-18.0 


16:56:48 


-40:14:26 


16533-4009 


7.8 


4.1 


NGG 6334F 


-10.4 " 


17:20:53 


-35:47:01 


17175-3544 


4.4 


5.5 


NGG 6334E 


e 


17:20:53 


-35:47:01 


17175-3544 


4.4 


5.5 



Table 1. Observed Sources, a) Phillips et al. 119981 : b) Walsh et al. 119981 : c) EUingsen et al. 120031 : d) Ellingsen ct al. 119961 . e) 

Observed within the same primary beam as NGC 6334F. Tabulated right ascension and declination values refer to the pointing centre 
(usually the methanol maser site), which is usually close to, but not necessarily coincident with the HII region centre. 




Figure 1. Gaussian least-squares fits to recombination line data. Error bars indicate the average noise level in each channel. Detection 
for source G309.92+0.48 is marginal, and a best-fit Gaussian is shown only for illustrative purposes. 



4 Shabala, Ellingsen, Kurtz & Forbes 




ingulnr Offset Right Aicensioii [J2000) 

Figure 2. H91q moment maps for G318.95— 0.20. Position-velocity plots are shown in part a (left plot). Contours are at 20 to 90% of 
maximum emission of 45.1 mjy beam~^ spaced by 10%. The lowest contour of 9.0 mjy beam~^ corresponds to 2.9(t. Velocity distribution 
across the region is shown in grey scale in part b (right plot). Also plotted in part b are continuum contours at 10 to 90% of maximum 
emission of 0.724 Jy beam-l spaced by 10%. From the 90° cross-cut (corresponding to a cut in RA), the slight western extension seems 
to move at a different characteristic velocity and is therefore probably not associated with other components, or is perhaps related to 
an outflow. A south-eastern extension is consistent with observed elongation in that direction on the continuum plot, thus the compact 
emission and an extension to the south-east are likely to be associated. Velocity slices suggest the amount of extended emission around 
this source is larger than for the others. The continuum and RRL emission are 2' distant from the methanol maser, suggesting that this 
region may be older than the other sources in our sample. 
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Figure 3. H91a moment maps for G328. 81-1-0. 63. The plots are as in Figure|2] Position-velocity contours are at 25 to 85% of maximum 
emission of 64.1 mjy beam~^ spaced by 10%, with the lowest contour of 16.0 mjy beam~^ corresponding to 4.1(t. Continuum contours 
are at 10 to 90% of peak emission of 1.46 Jy beam~^ spaced by 10%. Two distinct components are seen to the north and south of 
the compact core. Both are found to move at velocities close to those of the core, and are connected in the velocity slices. Hence the 
extended emission in this source is likely to be associated with the compact component. The north component coincides with the peak 
continuum emission observed in the continuum map, while the south extension corresponds to the cometary tail (see Figure 4 of ESK05). 
An additional north-west component separate from the rest of the emission also appears at approximately -35 km s~^. 
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Source 


^peak 


Sl 


FWHM 


Sc 






Tc 




source size 




(km s-l) 


(m.Ty) 


(km s-i) 


(.ly) 


(arcsec) 


(arcsec) 


(K) 


(K) 


(arcsec) 


G 308.92+0.12 


-44.4 ± 0.8 


19.9 ± 1.3 


34.7 ± 4.6 


0.263 ± 0.0002 


7.4 X 6.1 


7.3 X 6.8 


3900 ± 20 


8200 ± 1400 


8.6 X 4.0 


G 309.92+0.48 








0.670 ± 0.0003 


7.6 X 6.3 


7.4 X 7.0 






8.4 X 4.0 


G 318.95-0.20 


-29.0 ± 0.4 


45.1 ± 0.7 


26.4 ± 1.4 


0.724 ± 0.001 


8.1 X 6.4 


8.1 X 6.8 


9630 ± 40 


12600 ± 800 


7.7 X 4.0 


G 328.81+0.63 


-44.5 ± 0.3 


64.1 ± 1.2 


33.2 ± 0.9 


1.47 ± 0.0005 


8.2 X 6.1 


8.5 X 7.0 


18370 ± 80 


12900 ± 500 


6.6 X 4.0 


NGC 6334F 


-5.3 ± 0.2 


126.1 ± 1.6 


38.8 ± 0.8 


2.13 ± 0.06 


11.7 X 7.4 


11.3 X 7.6 


18300 ± 600 


10100 ± 500 


5.0 X 4.0 


NGG 6334E 


-4.0 ± 0.2 


57.1 ± 2.1 


23.4 ± 1.3 


0.417 ± 0.002 


11.7 X 7.4 


11.3 X 7.6 


3580 ± 30 


7500 ± 600 


5.0 X 4.0 



Table 2. H91a non-linear least squares Gaussian fits and derived parameters. Source sizes are FWZI. Sc values are taken from ESK05. 
Errors for the fits are estimated from Monte Carlo simulations (see text). Ypeak^ Sl, and FWHM are obtained from the Gaussian fits. 
The electron temperature, Te and continuum brightness temperature, Tc are derived as outlined in Section l3.2.1l For G 309.92+0.48 we 
have only a marginal detection, and hence no formal fit; as a consequence no is listed and no Te estimate can be made. 



PA SO degrees NCC: 6334 E 
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Figure 4. H91« moment maps for NGC 6334 E. The plots are as in Figurc|21 Position-velocity contours are at 35 to 95% of maximum 
emission of 57.3 mjy beam~^ spaced by 15%. The lowest contour of 20.1 mjy beam~^ corresponds to 3.6cr. Continuum contours are at 
82 and 90% of peak emission. The extended and central components of the RRL emission, separated by about 20 arcseconds on the sky, 
are practically coincident in the 60° cut shown. All cross-cuts suggest these move at a very similar velocity, and it is thus likely that they 
are associated. The level of background continuum emission is rather high in this source, resulting in poor image and p-v plot quality 
compared with other sources. High sensitivity VLA observations show NGC 6334 E to have a shell-like morphology with a weak, central 
compact source ICarral et al. 2002J . This structure is not well-imaged in the ATCA observations presented here. 



emission is observed around G 308.92+0.12, and almost none 
at all around G 309.92+0.48. These results are consistent 
with the continuum observations of ESK05, and in all cases 
suggest association between the compact and more extended 
components. 



3.2 Derived Source Parameters 

3.2.1 Stellar ionising flux and electron temperatures 

Assuming local thermodynamic equilibrium (LTE) , the elec- 
tron temperature Te is related to the continuum and recom- 
bination line brightness temperatures Tc and Tl, and the 
line Full Width at Half Maximum (FWHM) AV in km s'S 
by 



Although unlikely, line-of-sight effects cannot be ruled 
out, as illustrated by the star-forming complex NGC 6334. 
The E and F components of this complex are well-known 
to be separate star-forming regions, separated by more than 
one arcminute on the sky. However, their systemic velocities 
are very similar, with the best fit to the E component RRL 
giving a peak velocity of —4.0 ± 1.0 km s~^, compared with 
—5.3 ±0.6 km s~^ for the F component. For this reason, the 
two regions would have been extremely difficult to distin- 
guish had they been superimposed along our line of sight, 
rather than clearly separated on the sky. This suggests that 
all moment analysis results should be treated with a degree 
of caution. 



TlAV 

Tc 



6983u' 



where y*" is the fractional abundance of He"*" by number 
(McGcc fc Newton 198H . Assuming Y+ = 0.1 and rear- 
ranging yields 



6348 i/^'^Tc 
TlAV 



(1) 



The continuum and recombination line brightness tem- 
peratures are related to the peak flux density S (in Jy) and 
beam solid angle fl via 
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FA GQ degrees 



NGC 6334 F 




Angular Offset 



?*30"5e= S5" 54= 53' 

Right Ascenston U2000) 



Figure 5. RRL H 91 a moment maps for NGC 6334 F. The plots are as in Figure |2] Position-velocity contours are at 20 to 90% 
of maximum emission of 126.1 mjy beam~^ spaced by 10%. The lowest contour of 25.2 mjy beam~^ corresponds to 4.5f7. Continuum 
contours are at 10 to 90% of maximum emission of 2.13 Jy beam~^ spaced by 10%. Arcsecond-scale emission around the core is seen to the 
east and south of the emission centre. The southern component appears at velocities consistent with the -5.33 km s~^ systemic velocity 
of the region and represents the tail of the cometary region. The eastern component, however, peaks at velocities of about -f20 km s~^, 
which is significantly different from systemic region velocity and the maser velocities of approximately -10 km s~^. Therefore, it does 
not appear to be associated with the compact region. This result is corroborated by the lack of significant eastern extension in the 
superimposed continuum contours, or the high-resolution continuum maps of ESK05. This emission is in fact known to be coincident 
with a massive, poorly-coUimated bipolar outflow IBachiller &: Cernicharo 19901 Dackson et al. 19881 rHe" Free et aL 1995bl . 



[ Bright — 



S X 10-^'^ 



(2) 



The continuum observations were made at uc = 8.64 GHz, 
while the H91a rest frequency is i/l = 8.584 GHz. Te 
can then be determined by combining equations Q and |21 
In general, optical depth, pressure broadening, and stim- 
ulated emission must be accounted for in RRL analysis 
lUtoelfsema fc Goss 1992t . These non-LTE effects can be sig- 
nificant for the H91a line. For expected emission measures 
of ~ 10^ pc cm~® we use a correction factor of ~ 1.3 
llShaver et ai. 1983ll . The values of Te derived above must 
be scaled by this factor to account for non-LTE effects. The 
resulting values for Te are shown in Table |5| They are some- 
what higher than values obtained with single-dish obser- 
vations l |Gasweii fc Haynes 1987^ . The Te uncertainties for 
most sources are quite large due to uncertainties in the H91a 
fit parameters. 

We also calculate the ionising flux and stellar spectral 
type for each source by using the observed continuum fluxes 
of ESK05 and derived electron temperatures via the stan- 
dard Schraml & Mezger Ill9b9t argument. These values are 
given in Table El together with the corresponding spectral 
types from Panagia 1197311 . 



3.2.2 Emission Measures 

The continuum brightness temperature Tc is related to elec- 
tron temperature Te via the optical depth t, 

Tc ^Te{l- e"^) 

Rearranging this expression, the optical depth is given by 



V Te 



From equation |21 Tc and therefore r can be evaluated for 
each source using the derived Sc values and beam solid an- 
gles given in Table |5| From this, the peak emission measure 
for each source can be derived using 



EAT 



peak — 



8.235 X lQ-^a{v,Te)Te 



(4) 



(3) 



Here, a{v,Te) is a correction factor of order one adopted 
from Mezger & Henderson 119671 . Values of r and peak 
emission measure calculated for each source are shown in 
Table El 

The peak emission measures thus derived were used to 
estimate the angular and physical sizes of the ultra-compact 
components of the regions. Taking ultra-compact regions 
to have emission measures in excess of 10^ pc cm~®, the 
cutoff emission measure fraction was determined for each 
source. This was defined as the lowest contour level for 
which emission measure exceeds 10^ pc cm~'', and is given 
by cutoff— lO"^ / EMpeak. The location of the closest contour 
in the 6 km continuum images of ESK05 then determined 
the size of the observed ultra-compact region, given in Ta- 
ble El For the purposes of comparison with our models, the 
sources G 308.91+0.12, G 309.92+0.48 and G 318.95-0.20 
were considered to be spherical; while cometary sources 
G 328.81+0.63 and NGC6334F were modeled with the star 
offset from the centre of the spherical density distribution. 



4 MODELING 

Detailed numerical modeling will certainly be required to 
address the nature of extended emission associated with 
UC Hll regions. In this section, we present a simple, semi- 
quantitative model. Franco et al. 12000al I2000bt and Kim 
& Koo (IJOOl) suggest that the ambient density struc- 
ture is the primary factor determining Hll region sizes 
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Table 3. Ultra-compact component sizes derived from calculated emission measures and continuum contours of ESK05, and kinematic 
distances to the sources. Also given are the ionising flux for each source, derived from electron temperatures and continuum flux densities 
jSchr aml &; Mezgcr 1969 I of Table |3 and the corresponding spectral type (Panag ia 1973t . The distances to sources are from a) Phillips 
et al. (,1998 I; b) Caswell & Hayncs |1987|; a^nd EUingsen et al. (iI996l. Optical depths and peak emission measures are also given 
where appropriate. In the case of source G 309.92+0.48, Caswell & Haynes (19871 electron temperature scaled for non-LTE effects is 
used to derive EMp^^j^. A scaled Caswell & Haynes 119871 electron temperature is also quoted for source G 339.88 — 1.26. Continuum 
temperatures for sources G 328.81+0.63 and NGC 6334F are comparable with their electron temperatures, and we therefore adopt a 
lower limit of r = 1 for these sources. Correction factors a are interpolated from Mezger &; Henderson 119671 values, and range from 
0.9762 to 0.9940. The cutoff fractions have been determined by comparing the derived peak emission measures with the canonical UC Hll 
region value of 10^ pc cm~^. Large uncertainties associated with calculated Te and EMp^^k values arise due to propagation of smaller 
uncertainties through calculations. Peak emission measures predicted by the model outlined in Section |1] are also given; these should be 
treated as order-of-magnitude estimates. 



Source 


Spectral 


log5. 


Kin. Dist. 


Te 


T 


EMpeafe/10^ 


Cutoff 


UC size 


UC size 


predicted EMp^„fc/10^ 




type 


(s-i) 


(kpc) 


(K) 




(pc cm~^) 




(arcsec) 


(pc) 


(pc cm^^) 


G308.92+0.12 


BO 


47.01 


5.2" 


8200 ± 1400 


0.65 ± 0.12 


41 + 17 


0.024 


9 


0.21 


3.1 


G309.92+0.48 


7.5 


48.52 


5.3" 


12200 


1.54 + 0.40 


56 + 34 


0.018 


3 


0.078 


5.3 


G318.95-0.20 


BO 


47.72 


2.0'' 


12600 ± 800 


1.44 + 0.10 


56 + 8 


0.018 


13 


0.126 


2.6 


G328.81+0.63 


08 


48.39 


3.0'' 


12900 ± 500 


> 1 


41 + 4 


0.024 


6 X 11 


0.087 X 0.156 


4.3 


G336.40-0.25 


B0.5 


46.74 


5.2" 


4800 














G339.88-1.26 


B0.5 


45.99 


3.0^= 


10000 














G345.01+1.79 


BO 


47.11 


1.7^^ 


10000 














NGC 6334P 


09 


48.00 


1.7^= 


10100 ± 500 


> 1 


30 + 7 


0.034 


5x8 


0.041 X 0.066 


4.4 


NGC 6334E 


BO^ 


47.13 


1.7'= 


7500 ± 600 


0.65 + 0.07 


13 + 3 


0.079 






3.3 



and morphologies, thereby implying the need for more re- 
alistic ambient density representation; e.g. Franco et al. 
11199m . However, in the present work our focus is on the 
apparent association between the ultra-compact compo- 
nents and more diffuse arcsecond-scale extended emission 
IWood fc Churchwell 1989llKurtz et al. 1994L and we show 
that this can be explained in an order of magnitude argu- 
ment by a hierarchical density model. 

The density structure of star-forming cores is an 
important modeling parameter. Numerous studies have 
shown that in low-mass star-forming clouds the density 
structure on large scales (> 1 pc) is well- fit by power- 
law distributions n oc r^^. In high-mass star forma- 
tion regions the exponent of the density power-law flat- 
tens significantly for more evolved objects, such as Hll 
regions llBeuther et al. 20021 IHatcheU fc van der Tak 20031 
Ivan der Tak et al. 200011 . High-mass cores are less well- fit by 
single power-laws, and show a tendency towards clumpy sub- 
structure, possibly with the clumps embedded within over- 
all gradients (.Beuther et al. 2002. Evans et al. 1999j . Given 
the observational uncertainty regarding the magnitude of 
the density gradients and the scales over which they apply, 
we have ignored them in our modeling in favour of a simple 
approximation of a series of concentric spherical gas clumps. 

We assume that the star forms within a hot core {R = 
0.1 pc, To = 200 K, no = 10^ cm"''), that is located within a 
molecular clump of radius R — 0.3 pc, having molecular gas 
temperature To = 50 K and density no = 10^ cm~^, which 
itself lies within a molecular cloud {R > 0.3 pc, To = 25 K, 
no = lO' cm~^). The physical characteristics for the inte- 
rior hot core are taken from Churchwell (|5002), while those 
for the intermediate molecular clump are given by Cesaroni 
et al. II1991II and Garay & Lizano (1999 J, and for the exterior 
molecular cloud we used the parameters given by Church- 
well 119991 . We note that the values we use to define hot 
cores and molecular clumps are indicative only and differ 
slightly from those used by Kim & Koo 1200111 . 



In the simple model of Hll region evolution the radius of 
the expanding region is given as a function of time in terms 
of the initial Stromgren radius Rs and sound speed in the 
ionised gas by l |Dyson fc Williams "1980^ 

dR{t) _ ^ f Rit)\-'^' (5) 



dt 



The Stromgren radius is given as 

R.= i^^"''' 
\47rng/?2 

and assuming the strong-shock limit for the expansion 
following the (instantaneous) formation of the Stromgren 
sphere, we have 



Rit) = R,(^1 + 



ARs 



4.1 Model characteristics 

^.1.1 Thermal and turbulent pressure 

Given the higher thermal pressures that we now know to 
exist in molecular cores, Hll regions produced by 09 or later 
stars may still be ultra-compact when they reach pressure 
equilibrium with their surroundings IDe Free et al. 1995all . 
The non-thermal broadening of molecular lines in high-mass 
star forming regions suggests that turbulence is present, with 
velocities of the order of 2 km s~^ IXie et al. 19961 . The 
resulting additional turbulent pressure pturb = noniH^Vturb^ 
given in terms of the molecular hydrogen mass and the 
turbulent velocity Vturb in the surrounding medium, may 
act to restrict Hll region expansion. 

Using equation |K| we can compare the relative contri- 
butions of the expanding ionisation front (I-F) and tur- 
bulence in the ambient medium to the energy balance. 



/ dR \' 



-3/2 
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The sound speed in the ionised gas is = y ^^^^ ~ 

12.9 km s~^ for an electron temperature of 10 000 K. Taking 
an initial Stromgren radius of 0.02 pc, and UC region radius 
of 0.1 pc, we have R ^ bRs, and ^ 3.7. Thus, pho- 

toionisation energy nominally dominates (for Rk. Rs) but is 
of the same order as the turbulent energy. Turbulent veloci- 
ties greater than 2 km s~^ could shift the balance in favour 
of turbulence. Moreover, as the expansion proceeds, the I-F 
energy dominance will die off, as R grows well beyond Rs- 

4.1.2 Density structure in ionised regions 

Low-density extended emission on arcminute scales is 
observed near many UC Ifll regions (Kurt z et al. 1999 i 
IKim fc Koo 20011 [Ehingse n, Shabala & Kurt z 20051 . By 
comparison, as shown in Section 13.1.11 we observe emis- 
sion on arcsecond scales around the UC cores, consistent 
with other observations (e.g. Wood & Churchwell 1989; 
Kurtz et al. 1994). Inhomogeneous ambient density struc- 
ture can explain this IILi et al. 2005t . Non-uniformity within 
the ionised region can also arise if the expansion velocity of 
the ionisation and shock fronts is much greater than the 
sound speed, a condition that occurs early in the Hll region 
expansion phase. 

The expansion velocity of an Hll region slows with time, 
and is of the order of the sound speed when the region 
reaches pressure equilibrium. The diffusion timescale as the 
region expands into a molecular clump is of the order of the 
sound-crossing time t ~ . Taking ump 

~ 0.15 pc 

and Oi ~ 12.9 km s~^ as before gives t ~ 1.5 x 10* years. 
This is a significant fraction of an UC Hll region lifetime 
of ~ 10^ years, and hence the ionised gas density cannot be 
considered uniform in all cases. This situation is further am- 
plified by the presence of density inhomogeneities. Clearly, 
to model Hll regions properly, a full hydrodynamical treat- 
ment of the problem is required. Such modeling is beyond 
the scope of this paper, which purports only to offer a semi- 
quantitative plausibility argument. 

4.2 Comparison with observations 

Apart from NGC 6334E which happened to be in the same 
field of view as NGC 6334F, the nine regions presented here 
were selected for the presence of 6.7-GHz methanol maser 
emission. These masers are thought to correspond to a rel- 
atively short evolutionary phase that ends soon after the 
formation of the UCHll region (see ESK05). Recombination 
line analysis of extended emission around the majority of our 
sample shows that it is associated with the compact emission 
and thus the two must be considered together. We have com- 
pared the predictions from our model (compiled in Table l^l 
with the data for spherical (G 308.92-1-0.12, G 309.92-f0.48 
and G 318.95-0.20) and cometary (G 328.81-1-0.63) sources. 
The cometary source G 328.81-1-0.63 was modeled by posi- 
tioning the ionising star 0.096 pc from the hot core centre. 
In all cases, observed region sizes agree within a factor of 
a few with predicted pressure equilibrium values. However, 
the predicted peak emission measures are consistently more 
than an order of magnitude less than the observed values 
(see Table 1^. This discrepancy can be explained by the 



presence of significant amounts of ionised gas around the 
UC region on arcsecond scales, consistent with observational 
results of Section [3. 1.11 and discussed in more detail below. 
The remaining sources in our sample, particularly those with 
complex morphologies, will require more detailed modeling 
than is considered here. 

For a spherical Hll region the distance of a site line 
from the centre of the region is d = rcosO. The emis- 
sion measure at this distance is obtained by traversing a 
length 2rsin9. For uniform electron density, we then have 
EM{d)/EMp^ak = {2nlrsine) / (2n^r) = sinO. Thus 

d f f EM{d) y 

The resultant theoretical contours can then be com- 
pared with observations. The ultra-compact components of 
sources G 308. 92-^0.12, G 309. 92-^0.48 and G 318.95-0.20 
were largely unresolved in the 750D array results presented 
in ESK05, and high resolution images made with the 6 km 
array (also presented in ESK05) were used instead. The 
synthesized beam FWHM was taken as 1.2 arcseconds for 
all three sources. For each source, the resulting beam was 
then convolved with theoretical contours using the Table 
1^ ultra-compact region sizes. Figure |n] shows the theoreti- 
cal map thus obtained, together with a high resolution im- 
age, for G308.92-f0.12. Table gl shows the theoretical con- 
tour diameter and observed major and minor axes for each 
contour of each source. The final three columns give the 
observed/theoretical ratios for the two axes, as well as a 
geometric mean of the two for sources G 308.92-1-0.12 and 
G 309. 924-0.48. The non-spherical nature of G 318.95-0.20 
(much more so than the other two sources) means that we 
have only given the major axis values for this source. 

In all three sources the observed/theoretical ratios de- 
crease as we approach peak emission. Thus the observed 
contours are slightly denser near source core, implying 
the presence of density gradients in the ionised gas. For 
G 309.92-1-0.48 and especially G 308.92-1-0.12 the ratios are 
close to constant, suggesting almost uniform region density 
and thus that these sources are close to pressure equilib- 
rium. This is again consistent with a lack of lower density 
gas observed around their UC cores. By comparison, the 
large spread of observed/theoretical ratios in G 318.95—0.20 
suggests a much steeper ionised gas density gradient in this 
source, in keeping with observations of significant arcsecond- 
scale emission around its UC core. Accounting for this den- 
sity gradient would raise the predicted peak emission mea- 
sure and thus address the discrepancy between model and 
observations discussed in the previous section. 

The ratios given in Table |1] are less than one for all 
three sources, indicating that the observed Hll region sizes 
are smaller than model predictions. This could be due to 
the Hll regions not being in pressure equilibrium with the 
ambient medium — an idea consistent with their young ages 
deduced from ma ser observations (lEllingsen et al. 1996| 
[Phillips et al. 1998|IT)e Buizer et al. 200 2^ and also the fact 
that this ratio is closer to one for sources G 308.92-1-0.12 
and G 309.92-1-0.48 which exhibit a more uniform density 
structure. G 309.92-1-0.48 is unresolved in the 750D array, 
and this is likely the main reason for the departure of ob- 
served contours from model predictions. Overestimates of 
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Figure 6. Observed (left) and theoretical (right) emission measure contours for G 308.92+0.12. Both sets of contours are at 32, 47, 62, 
77 and 92% of the peak emission measure. The box plotted together with the observed contours encloses the lowest (32%) theoretical 
contour. The predicted and observed contours are of similar physical size, with some evidence for compact core- extended arcsecond-scale 
emission structure. 



Table 4. Theoretical and observed sizes for the three sources modeled as spherical. The observed/theoretical major and minor axis ratios 
are given for each contour of each source. The geometric means of the two are also given for sources G 308.92+0.12 and G 309.92+0.48. 
In source G 308.92+0.12, the major axis ratio is almost constant and close to one for most contours, indicating good agreement between 
model and observations. By comparison, the major axis ratios are less and depart significantly from unity in source G 318.95— 0.20, 
consistent with extended emission observed around this source (see text). The contours for source G 309.92+0.48, although in reasonable 
agreement with the model, are most likely unresolved. The decrease in the observed/theoretical ratio at higher emission measure contours 
provides support for existence of ionised density gradients in the UC cores of all three Hll regions. 



Source 


EM contour 


Theoretical diameter 


Major 


Minor 


Major axis 


Minor axis 


Geometric 




(%) 


(arcsec) 


axis (arcsec) 


axis (arcsec) 


ratio 


ratio 


mean ratio 


G 308.92+0.12 


32 


4.36 


4.28 


2.69 


0.98 


0.62 


0.78 




47 


3.80 


3.61 


1.89 


0.95 


0.47 


0.67 




62 


3.29 


3.11 


1.42 


0.95 


0.43 


0.64 




77 


2.65 


2.49 


0.94 


0.94 


0.35 


0.57 




92 


1.66 


1.49 


0.52 


0.90 


0.31 


0.53 


G 309.92+0.48 


11 




2.85 


2.34 










31 




2.05 


1.72 










51 


2.22 


1.77 


1.56 


0.80 


0.70 


0.75 




71 


1.67 


1.31 


1.14 


0.78 


0.68 


0.73 




91 


0.93 


0.56 


0.47 


0.60 


0.51 


0.55 


G 318.95-0.20 


32.5 


12.58 


8.22 




0.65 








52.5 


11.31 


6.22 




0.55 








72.5 


9.23 


3.24 




0.35 








92.5 


5.08 


1.25 




0.25 







stellar spectral types are another possible reason for the ob- 
served/predicted ratios being less than one, although this 
is less likely as radio observations typically underestimate 
spectral types due to dust absorption. Other confinement 
mechanisms may also play a role. Evidence for non-thermal 
broadening in the Gaussian profiles of Figure^lends further 
support to this scenario. 

The above analysis is applicable to optically thin Hll 
regions. If we instead had a constant continuum brightness 
temperature (as would be expected for an optically thick 
source), the theoretical emission measures would be more 



uniform around the source core, providing an even greater 
discrepancy between predictions and observations. 

5 DISCUSSION 

5.1 Emission Measures 

The fact that some fraction of ionising photons is absorbed 
by dust suggests that our measured peak emission mea- 
sures, which are already too high to be explained by con- 
stant density models, are underestimates. This effect can 
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largely be ignored however, as the attenuation factor is 
(1 — fY^^ , where / is the fraction of photons absorbed by 
dust ({Franco et al. 1990|l . which for / ~ 0.9 results in a de- 
crease in the predicted peak emission measure by only a 
factor of two. 

The large uncertainties associated with the electron 
temperatures we derived affect the calculated peak emis- 
sion measures both directly and through the value of optical 
depth r in equation The logarithmic dependence of r on 
Te will have a greater effect on the estimated peak emis- 
sion measure value than directly through the power-law Te 
dependence in the denominator of equation 2] Hence the 
higher electron temperature estimates we have calculated 
compared to Caswell & Haynes 119871 suggest that, if any- 
thing, the estimated peak emission measures are likely to be 
underestimates. Typically, these effects largely cancel each 
other out, and in any case affect the derived peak emission 
measures by at most a factor of a few. 

5.2 Lifetimes 

The fundamental difference between pressure confinement 
and other models is that it predicts that in many cases the 
observed Hll regions are already in equilibrium with their 
surroundings, rather than still undergoing expansion. Our 
modeling suggests pressure equilibrium may be reached very 
quickly, with expansion taking place for only a fraction of the 
observed lifetimes of Hll regions. The lower Hll region age 
limits thus derived are given in Table|K| These are consistent 
with methanol maser emission being associated with very 
young massive stars. 

As discussed in section l4.1.1l turbulence may provide an 
additional confinement mechanism. Any non-isotropic na- 
ture of the turbulence (e.g., if it is magnetohydrodynamic; 
Garcia-Segura & Franco 1996) may also contribute to the 
non-spherical appearance of the resulting Hll region. Fur- 
ther investigation of this issue is warranted. 



6 CONCLUSIONS 

We have detected arcsecond-scale emission around UC Hll 
cores. Using region parameters derived from continuum and 
H91a recombination line data we show that although sim- 
ple models of expansion in hydrostatic equilibrium repro- 
duce the observed region sizes, their emission measures are 
significantly underestimated. This discrepancy can be ex- 
plained by the presence of density gradients in the ionised 
gas, consistent with young source ages and observations of 
the diffuse emission. 
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Table 5. Lower limits on ultra-compact component ages as predicted by the hierarchical ambient density structure model. These are 
obtained from the time required to reach pressure equilibrium for each source. 
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Minimum Age (years) 
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2 X 10* 


G 309.92+0.48 
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8 X 10* 
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